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ABSTRACT 

Double-barred galaxies are common in the local Universe, with approximately one 
third of barred spirals hosting an smaller, inner bar. Nested bars have been proposed 
as a mechanism for transporting gas to the very central regions of the galaxy, trigger 
star formation and contribute to the growth of the bulge. To test this idea, we per- 
form for the first time a detailed analysis of the photometry, kinematics and stellar 
populations of a double-barred galaxy: NGC 357. We find that this galaxy is either 
hosting a pseudobulge or a classical bulge together with an inner disc. We compare 
the relative mean luminosity-weighted age, metallicity and a-enhancement between 
the (pseudo)bulge, inner bar and outer bar, finding that the three structures are 
nearly coeval and old. Moreover, the bulge and inner bar present the same metallic- 
ity and overabundance, whereas the outer bar tends to be less metal-rich and more 
a-enhanced. These results point out that, rather than the classical secular scenario in 
which gas and star formation play a major role, the redistribution of the existing stars 
is driving the formation of the inner structures. 

Key words: galaxies: individual: NGC 357 — galaxies: photometry — galaxies: kine- 
matics and dynamics — galaxies: bulges — galaxies: structure — galaxies: evolution. 
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1 INTRODUCTION 

The structural evolution of disc galsixies and its timescale 
remain as open questions nowadays. In fact, it is still not 
clear if the formation of bulges is mostly driven by exter- 
nal processes, such as interactions or mergers, or by the 
internal secular evolution due to inflow of gas from the 
galactic disc to th e central regions, where it ma y trigger 
star formation (see iKormendv fc Kennicutdl20o3 . for a re- 
view). This second scenario implies a slow process that 
takes place when non-axysimmetric components (particu- 
larly bars) redistribute the angular momentum of the galaxy. 



* Based on observations carried out at the European Southern 
Observatory (ESO 70.B-0338). 
f adlcr@iac.es 

I vazdekis@iac.es 
§ jalfonso@iac.es 

If enricomaria.corsini@unipd.it 

II vpdebattista@gmail.com 



Given the large number of barred galaxies fo und i n the 
local Universe |Knapen et al. 20001: Eskridge et al, 200C; 



Marinova fc Jogeel l2007l : iMenendez-Delmestre et al.1 boOTI : 



Aguerri et al. 20091 ). and that the bar fraction remains sig- 



nificant till z~l (althoug h its evolution with redshift is still a 
matter of debate, see e.g. lElmegreen et al ]|2q04l : lJogee" et al.l 
l2004l : ISheth et al.ll200alCameron et al.ll2010h . it seems that 
internal secular processes might actually be an important 
driver of the evolution in spirals. If this is the case, the re- 
sulting kinematical and morphological changes could also 
lead to an evolution along the Hubble sequence from late r 
to earlier types, as sugg ested bv lFriedli fc Martinetl |l993). 

It has been theoretically demonstrated that the flow of 
gas through a single bar stops before reaching the galac- 
tic centre; however, a system of nested bars may efficiently 
transport the material to the nuclear regions, where it even 
would hel p to fuel nuclear activ i ty. Th is idea was first pro- 
posed by IShlosman et al. 1 1| 19891 . liggoD . and resulted in an 
increasing interest in galaxies hosting two or more bars. 
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Double-barred galaxies are common objects in the lo- 
cal Universe; at least one third of barred galaxies show an 
additional inner bar (Mgiscev 2001; Er win fc Spar ke 2002; 
iLaine et al.l |2002| ; lErwinll2004 ) and n ested bars have been 
found till z~0.15 ( Lisker et al.ll2006l ). Photometric studies 
have shown that t he two bars are randomly oriented rela- 
tive t o each other (|Friedli fc Martinetlll993l : IWozniak et alj 
Il995l ). suggesting that they rota te independently, a.s pre- 
dicted by numerica l simulations ijShlosman fc Helled [200^ : 
iDebattista fc ShenI [2OO7I ) and confirmed by observations 
i Corsini et al.lioosi ). Moreover, it seems there is no correla- 
tion between the presence of an inner bar and the Hub- 
ble type of the host galaxy or its main bar properties. 
lErwin fc Sparkd l|2002l ') find that inner bars extend to ~12% 
of the outer bar size, and that have in general smaller el- 
lipticities, but may be larger since the inner bar might be 
embedded in the light of the bulge. 

iMacieiewski fc Spark3 (|l99l[200ol ) study the dynamics 
of double-barred objects and develop a simple formalism to 
study the orbits of their stars based on the loop concept. 
Single bars are characterized by two fundamental frequen- 
cies (one for the bar and one for the free oscillations), allow- 
ing the existence of closed and periodic orbits. In fact, two 
orbital families stand out for a single-barred potential: the 
xi orbits, elongated along the bar major axis, and the X2 or- 
bits, perpendicular to it (|Contopoulos fc PapavannopoulosI 
Il980l ). A double-barred potential needs an additional fre- 
quency to be described, meaning that closed orbits are not 
possible. A loop is a family of orbits whose populating par- 
ticles return to the same curve when the two bars recover 
the relative position between them, although each individual 
particle is not in its initial position. This means that only 
the global shape of the orbit is preserved. 

The challenge of building long-lived double bars 
led to many numerical simulations, most of them get- 
ting galaxies where the outer bar appears first, the 
gas flows along it and then it is captured by the 
X2 orbits of the main bar, subsequently forming the 
inner bar. In this scenario the dissipative component 
plays a major role (|FYiedli fc Martinetl Il993l; Heller et ahl 
2OOII; IShlosman fc Helleij |2002| : iRautiainen et al.l |2002| ; 



Englmaier fc Shlosmanl |2004| ). However, some authors get 



double-barred systems without such a dissipational process. 
This is the case of iDebattista fc S hen (2007j), who perform 
coUisionless N-body simulations and obtain purely stellar in- 
ner bars. These inn er bars form before the apparition of the 
outer bars. Finally. iHeller et al.l (2007) generate for the first 
time a two bar system from a dark matter halo. In this case, 
both bars are gas-rich and transport material to the central 
regio ns, supporting the hypothesis of lShlosman et al] l| 19891 . 
I1990I ). 

With these pieces of evidence, it is clear that double- 
barred galaxies are the perfect benchmark to constrain the 
role of secular processes within galaxy nuclei. One way to 
proceed is to perform a detailed stellar population anal- 
ysis of these systems: if double bars are playing a ma- 
jor role in galaxy evolution, this fact will be refiected in 
the properties (star formation, age, and metallicity) of the 
galaxy components (bulge, inner bar, outer bar, and disc), 
resulting in a gradient of ages and metallicities between 
them. However, this is difficult to achieve due to the struc- 
tural complexity of galaxies with nested bars, as illustrated 



by the photometric properties mentioned above. To date, 
there are few stellar population studies involving double 
bars, all of them being projects devoted to the analysis of 
samples of barred galaxies in general, which include only 
some double-barred objects. For this reason, these works 
focus only on the main bar and they do not account for 
the inner bar. This is the case in the set of papers b y 
iPerez et al.l (l2007l.l2009l),|Perez fc Sanchez-Blazguezl (|201 j ). 
and ISanchez-Blazguez et al.l ( 2011^ ■ who study the stellar 
populations of bars and bulges. They find positive, negative, 
and even null metallicity gradients along the bars indepen- 
dently of the age distribution. The bulges hosted in barred 
galaxies tend to be more metal-rich than those in unbarred 
galaxies at a similar velocity dispersion. However, the impli- 
cations of these results within evolutionary scenarios is still 
unclear. 

Other works study double-barred objects from 
a pu rely k inematical point of view (|Emsellem et al.l 
200 ll; iMoiseev. Valdes. fc ChavushvanI |2004|: 



de Lorenzo-Caceres et al. 20081 ). de Lorenzo-Caceres et all 
1200^ present the 2D stellar velocity and velocity dispersion 
maps for a sample of four double-barred galaxies, based on 
observations with the SAURON integral field spectrograph. 
The analysis of high quality velocity dispersion maps 
revealed two local minima, located exactly at the ends 
of the inne r bar of each galaxy. By mean s of numerical 
simulations, Ide Lorenzo-Caceres et al] (|2008l ) conclude that 
these a-hollows appear because of the strong contrast 
between the velocity dispersions of the bulge and the inner 
bar. This result indicates that bars are dynamically cold 
components, with velocity dispersion values significantly 
lower than those of the hot, classical bulge. Moreover, the 
(j-hoUows are the kinematical signature of the presence of 
an inner bar. 

In this paper, we present the first detailed morphologi- 
cal, kinematical and stellar population analysis of the differ- 
ent structures that comprise a double-barred galaxy: bulge, 
inner bar and outer bar. We focus on the case of NGC 357. 
The work is organised as follows: in Section [2] we list the 
relevant properties of NGC 357 from the literature, while in 
Section|3]we summarise the observations and the data reduc- 
tion procedure. Then, we describe the morphology and kine- 
matical profiles in Sections 14.11 and 14.21 respectively. Section 
4.3l details the stellar population analysis, from the emission 
lines correction of the spectra to the detailed study of the 
bulge region and the comparison between structural compo- 
nents. A discussion of the implications of these results for 
the structure and possible formation scenarios of NGC 357 
can be found in Section [S] Finally, Section [S] sets out the 
main conclusions of this work. 



2 MAIN PROPERTIES OF NGC 357 



NGC 357 is classified a s SB(r)0/a bv[de Vaucouleurs et"aLl 
(1991") and as SBa by ^Sandage fc Tammannl (|l98ll ). This 
g alaxy appears in the catalogue of double-barred galaxies 
of lErwir] (|2004l ). who provides information about the prop- 
erties of the main and the inner bar of NGC 357, based on 
the analysis of infrared images from the Hubble Space Tele- 
scope. The two bars of NGC 357 are almost perpendicular, 
with the inner one almost parallel to the major axis of the 
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Table 1. Relevant properties of the structural components of 
NGC 357. 



NGC 357 




Inclination 


PA 


Semi-major axis 






(1) 


(2) 


(3) 


(4) 


Inner bar 




45° 


3.1 arcsec 


0.16 


Outer bar 




120° 


21 arcsec 


0.44 


Disc 


36° 


20° 







(1) galaxy inclination llAguerri et al.ll2005l ): 

(2) p osition angle (hereaf ter PA) of the ma j or axis of the galaxy 
(disc; lAguerri et"aLll2005^ and the two bars l lErwinll2003 '): 

(3) bar lengths, esti mated as th e semi-major axes of maximum 
isophotal ellipticity l|Erwinll2004h . and 

(4) maximum isophotal ellipticity of the bars {e = 1 — b/a), from 
lErwinI |20oi)- 

galaxy, as seen in Figure [T] We have specifically chosen this 
early-type spiral in order to avoid the presence of complex 
structural components and dust, which can make the stellar 
popul ation analys is very tricky. The galaxy distance is 31.6 
Mpc (|Erwinl l2004'). which corresponds to ~150 pc arcsec 
The relevant properties of NGC 357 are shown in Table [T] 

Although some studies claim this galaxy is isolated 
(e.g.. lGadotti et al.ll2007l ).l van den Berghl (|2002l ) finds it be- 
longs to a group with at least six other members. However, 
it shows no signatures of interaction with i ts closer com- 
panio ns. NGC 357 is classified as a LINER (|Gadotti et al.l 
l2007l ). 

The structural properties of NGC 357 have been 
studied photometric ally by several authors. In this way, 
lAguerri et al.l l|2005l ) perform a photometric decomposition 
of an /-band image of this galaxy along its outer bar major 
axis, getting the main photometric parameters of its struc- 
tural components, namely bulge, outer bar and disc; they 
do not take into account the inner bar since its contribution 
to the radial surface brightness along the main bar is neg- 
ligible. Their main conclusion is that the bulge of NGC 357 
follows the same fundamental plane than the ellipticals and 
other bulges of SO galaxies. 



3 OBSERVATIONS AND DATA REDUCTION 

The spectroscopic observations of NGC 357 were carried 
out with the 3.5-m New Technology Telescope (NTT) at 
the European Southern Observatory (ESO) in La Silla 
(Chile) on 5-13 October 2002. The ESO Multi-Mode Instru- 
ment (EMMI) was operated both in blue (BLMD) and red 
medium-dispersion spectroscopic (REMD) mode. 

The NTT mounted EMMI in BLMD using the grating 
No. 3 blazed at 3800 A with 1200 grooves mm ^ in first or- 
der in combination with a 1.3 arcsec x 5.5 arcmin slit. The 
detector was the No. 31 Tektronix TK1024 EB CCD with 
1024 X 1024 pixels of 24 x 24 fim'^. It yielded a wavelength 
coverage between about 3990 A and 4440 A with a reciprocal 
dispersion of 0.45 A pixel"^. The spatial scale was 0.37 arc- 
sec pixel"^. The instrumental resolution was 1.3 A (FWHM) 
corresponding to (Jinst ~ 40 km s~^. 

Six spectra of 45 min each were taken aligning the slit 
with the major axis of the inner bar (PA = 45°). Two more 



spectra of 45 min each were obtained aligning the slit with 
the major axis of the outer bar (PA — 120°). All the spectra 
were obtained using the guiding TV camera to center the 
slit on the galaxy nucleus. 

The NTT mounted EMMI in REMD using the grat- 
ing No. 6 blazed at 5200 A with 1200 grooves mm"'^ in first 
order. A 1.0 arcsec x 5.5 arcmin slit was adopted. The mo- 
saiced MIT/LL CCDs No. 62 and 63 with 2048 x 4096 pixels 
of 15 X 15 fj,m^ covered the wavelength range between about 
4800 A and 5460 A. The on-chip 2x2 pixel binning provided 
a reciprocal dispersion and spatial scale of 0.40 A pixeP^ 
and 0.332 arcsec pixel"^. The instrumental resolution was 
1.6 A (FWHM) corresponding to ainst ~ 30 km s"\ 

After centering the slit on the galaxy nucleus, two spec- 
tra of 30 min each were obtained along the disc major axis 
(PA = 20°). 

The range of the seeing FWHM during the observing 
runs was 0.6-1.4 arcsec as measured by the ESO Differen- 
tial Image Motion Monitor. A comparison lamp exposure 
was obtained after each object integration to allow accu- 
rate wavelength calibration. Quartz lamp and twilight sky 
flatfields were used to remove pixel-to-pixel variations and 
large-scale illumination patterns. Several G and K stars and 
spectrophotometric standard stars were observed with the 
same set-up to serve as templates in measuring the stellar 
kinematics and in fiux calibration, respectively. 

All the spectra were overscan and bias subtracted, flat- 
field corrected, corrected for bad pixels and columns, and 
wavelength calibrated using standard IRAfQ routines. The 
cosmic ray removal is a critical step since any residual might 
affect the spectral lines measured to derive the properties 
of the stellar popula tions. It was performed with the RE- 
DUCEME package (jCardiell [l999l ). that assures a careful 
and accurate inspection and interpolation of the spectra. We 
checked that the wavelength rebinning was done properly 
by computing the d ifference between the m easured and pre- 
dicted wavelengths (|Osterbrock et akll 19961 ) for the brightest 
night-sky emission lines in the observed spectral ranges. The 
resulting accuracy in the wavelength calibration is better 
than 5 km s~^. The spectra taken along the same axis were 
co-added using the centre of the stellar continuum as refer- 
ence. The contribution of the sky was determined from the 
outermost regions at the two edges of the resulting spectra, 
where the galaxy light was negligible, and then subtracted, 
giving a sky subtraction better than 1%. A one-dimensional 
sky-subtracted spectrum was obtained for each kinematical 
template star. The sky-subtracted spectra were flux cali- 
brated using the observed spectrophotometric stars as a ref- 
erence. Neither extinction nor dust corrections were applied 
for the blue spectra since the wavelength coverage is short 
and we are not interested in colour measurements. 

Finally, the galaxy spectra were binned along the spa- 
tial direction in order to assure a minimum signal-to-noise 
ratio (S/N hereafter) of ~ 20 A~^, sufficient to our kine- 
matic analysis. Such a minimum S/N is reached only in the 



^ Imaging Reduction and Analysis Facilities (IRAF) is dis- 
tributed by the National Optical Astronomy Observatories which 
are operated by the Association of Universities for Research in 
Astronomy (AURA) under cooperative agreement with the Na- 
tional Science Foundation. 
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Figure 1. r-band image of NGC 357 from the Sloan Digital Sky Survey l|York et al.l200d ). The left panel shows the whole galaxy whereas 
the right panel is a zoom of the inner ~20 arcsec. We have overplotted the isodensity contours and the directions of the major axis of 
the inner bar (green dot-dashed line), outer bar (green dashed line) and disc (white line). 



outermost radial bins, it increases at smaller radii reaching 
a maximum S/N~ 60 in the radial bin corresponding to 
the galajcy centre. The spatial binning procedure was care- 
fully carried out in order to assure the separation of the 
different structural components of NGC 357 (i.e., the bulge, 
inner, and outer bar). 



4 ANALYSIS AND RESULTS 
4.1 Photometry 

As a first approach to disentangling the morphology of 
NGC 357, we perform an ellipse fitting over a NICMOS2 
F160W image taken from the Hubble Legacy Archive (Prog. 
Id. 7330, PI John Mulchaey). For this purpose we use the 
IR AF task ELLI PSE, which follows the procedure outlined 
bv IJedrzeiewslH ((1987, ). Figure [2] shows the resulting pro- 
files of elUpticity, e, position angle, PA, and the fourth co- 
sine Fourier coefficient, a4. The presence of the bars causes 
two sharp changes in the ellipticity and PA parameters: one 
peak at ~3 arcsec related to the inner bar, and a disconti- 
nuity at ~25 arcsec due to the outer bar. These two values 
are in agreement with the bar lengths estimated bv lErwinI 
(|2004l ). The a4 profile also shows two peaks at the locations 
of the two bars. 04 is a measurement of the deviation of the 
isophotes from pure ellipses, so positive values correspond to 
discy isophotes wh e reas negative values to boxy isophotes 
l|Jedrzeiewskilll987l : [Kormendv et alll2009l ). Therefore, the 
two bars of NGC 357 present discy isophotes. 

We get a maximum isophotal ellipticity of e ~0.15 and 
e >0.4 for the inner and outer bars, respectively, and a 
PA~40° for the inner bar, PA~120° for the outer bar, and 
PA~15° for the disc. The se measureme nts a re also consis- 
tent w ith the estimates of lErwii] (|2004l ) and lAguerri et al] 
l|2005l ). 



4.2 Kinematics 

4-2.1 Measuring the line-of-sight velocity and velocity 
dispersion 

Two different techniques were used to derive the velocity and 
velocity dispersion profiles of NGC 357: a cross-correlation in 
the Fourier domain and a full spectrum fittin g in wavelength 
space . The cross-correlation procedure (e.g.. lTonrv fc David 
I1979I ) is as follows: among the observed velocity standards, 
we select those that are spectroscopically most similar to 
the galaxy. After correcting them to rest frame, we com- 
bine the stellar spectra to create an artificial velocity stan- 
dard star spectrum that will be used as a template for the 
cross-correlation. This standard spectrum is then broadened 
to several velocity dispersions and the results are corre- 
lated with the original star, obtaining for each case a cross- 
correlation peak with its full width half maximum (FWHM). 
Thus, we create a calibration table that provides the FWHM 
we expect if a spectrum with a given velocity dispersion is 
cross-correlated with our stellar template, so we can get the 
dispersion of the galaxy spectra by correlating them with 
the star. The line-of-sight velocity is a direct product of the 
correlation since the template is at rest. 

For the full spectrum fitting we have made use of the pe- 
nalized pixel fitting (he r eafter pPXF) routine developed by 
ICappellari fc Emselleiiil (|2004l ). The spectral regions which 
are potentially affected by emission lines are previously 
masked on each galaxy absorption spectrum, which is then 
fitted with a linear combination of a well-selected subsample 
of the single stellar | 3opulation models (hereafter SSPs) of 
'Vazdckis ct al. (zOlO), previously convolved with a line-of- 
sight velocity distribution (hereafter LOSVD). This LOSVD 
is parametrized by an expansion in Gauss-He rmite functions 
(|Gerhardlll993l :l van der Marel fc Franxlll993l ): since the S/N 
requirements for this procedure get more demanding as the 
order of the desired LOSVD moments increases, we decided 
to focus only on the velocity and velocity dispersion to avoid 
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Figure 2. EUipticity (left panel), position angle (middle panel), and fourth cosine Fourier coefficient (right panel) profiles that result 
from fitting ellipses to the isophotes of NGC 357 in the NICMOS2 F160W image. The semi-major axes of the fitted ellipses are given in 
logarithmic scale. The vertical orange dot-dashed and blue dashed lines indicate the lengths of the inner and outer bars, respectively. 



the loss of spatial resolution caused by a wider radial bin- 
ning. 

The two methods used for deriving the kinematics of 
NGC 357 are completely independent and even work in 
different parameter spaces. In order to enhance the dif- 
ferences between both techniques, we explore two differ- 
ent ways that offer complementary advantages: whereas 
the SSPs help to alleviate t he template mismatch problem 
l|Falc6n-Barroso et al.ll2003l ). the artificial template was ob- 
tained by stellar spectra, which were observed and reduced 
exactly in the same way as the galaxy spectra, thus mini- 
mizing the possibility of introducing instrumental effects. 

The results obtained from both methods are fully con- 
sistent, with differences smaller than the corresponding error 
bars. These differences range from few kms"'^ in the central 
regions to a maximum of 10 kms^^ in the farthest observed 
radii. Figure |3] shows the final line-of-sight velocity (after 
subtracting the systemic velocity) and velocity dispersion 
profiles obtained with pPXF along the major axis of the 
inner bar, outer bar, and disc. 

4-2.2 Revealing the inner and outer bars 

The presence of the inner bar is clearly unveiled by the shape 
of the rotation curve (Figure [31 left panels), that follows the 
same characteristic S- or double-hump profile pr eviously ob- 
served for inner bars by lEmsellem et all (1200 ll ). This pro- 
file is characterized by an inner steep slope so the rotation 
reaches a local maximum, it slightly drops to a local mini- 
mum and then it raises again slowly further out. Since the 
two bars of NGC 357 are almost perpendicular, with the 
inner one almost parallel to the major axis of the galaxy, 
we can conclude that the double-hump profile in Figure [3] is 
completely due to the presence of the inne r bar. This kind of 
rotati on cur ve was theoretically derive d bv | Chung fc Bureaul 
l|2004) and [Bureau fc Athanassoulal (|2005l ) by means of 
N-body simulations of single-barred galaxies; in this last 
work, the authors show that the double-hump structure is 
more enhanced in galaxies with strong bars seen end-on, 
whereEis the effect is weakened for larger viewing angles 
or weaker bars. Although the double-hump rotation curve 
is clearly visible in the case of the inner bar of NGC 357, 



suggesting the bar is strong, the maximum isophotal ellip- 
ticity is only emax ~0.15, which indicates that it is actu- 
ally a very weak bar. Indeed, the maximum ellipticity of 
the isoph otes in a bar region co rrelates well with the bar 
strength (jLaurikainen et al1l2002l ). 

The velocity dispersion profile along the inner bar shows 
three local mini ma: a central dr o p cor responding to the so- 
called (J -drop of lEmsellem et ahl (|200ll ). and two symmetric 
decreases exactly at the edges of the bar (±4 arcsec) that 
can b e identified as the a-hoUows l|de Lorenzo-Caceres et ahl 
l2008l ). known to be signatures of an inner bar. These hol- 
lows have an amplitude of ~20 kms~^, co nsistent with the 
findings of lde Lorenzo-Caceres et al.l l|2008l ). who found am- 
plitudes ranging from 10 to 40 kms~^. Further discussion on 
these features and their different origins is found in Section 
[5] As expected, the velocity dispersion decreases outwards. 

The outer bar of NGC 357 (Figure [S] middle panels) is 
almost perpendicular to the major axis of the galaxy, which 
we assume to be the kinematical line of nodes. This means 
that the line-of-sight velocity values along the outer bar are 
expected to be very low, as it is shown in Figure O However, 
a clear kinematical decoupled profile appears in the central 
region. This structure rotates faster than its surroundings 
and extends ~ ±2 arcsec, which corresponds to the radial re- 
gion of the cr-drop observed in the velocity dispersion profile. 
A closer inspection to the inner bar velocity profile reveals 
the same decoupling, since two steep humps appear inside 
the double-hump profile, exactly at ~ ±2 arcsec. Unfortu- 
nately, the low S/N level of the spectra along the direction 
of the outer bar of NGC 357 only allows us to analyse the 
kinematics in the inner ~ ±10 arcsec, so we cannot explore 
the whole main bar to see if other signatures, as for example 
the (j-hoUows or the double-hump profile, are present. 

The kinematics along the major axis of the disc is shown 
in Figure [3] (right panels) . In this case, we find as expected 
the maximum line-of-sight velocity values and the same cen- 
tral signatures found for the other directions: the decoupling 
in the velocity profile and the cr-drop. In Section [S] we dis- 
cuss the possible structures that might be causing these sig- 
natures. 
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Figure 3. Linc-of-sight velocity (after tiie systemic velocity subtraction; top panels) and velocity dispersion (bottom panels) profiles as 
a function of the radius for the spectra obtained along the major axis of the inner bar (left panels), outer bar (middle panels) and disc 
(right panels) of NGC 357. 



4.3 Stellar populations 

The aim of this stellar population analysis is to recover the 
relative ages, metallicities, and formation timescales of the 
different structural components of NGC 357: bulge, inner 
bar, and outer bar. Unfortunately, the quality of the data, 
although very high, still is not sufficient to analyse the disc, 
which requires spectra taken with 10 m-class telescopes. To 
perform this detail study, we have summed up the individual 
spectra, previously corrected for their kinematics (velocity 
and velocity dispersion), to reach the desired S/N within 
each structural component. For this purpose, we only add 
those spectra corresponding to the regions where a given 
component is clearly dominating the total luminosity: the 
inner ±1 arcsec for the bulge, from ±2 arcsec to ±4 arc- 
sec for the inner bar, and from ±4 arcsec to ~ ±15 arcsec 
(limited in this case by S /N requirements) for the outer bar. 
Since the spectra along the major axes of the two bars and 
the disc always cover the central region, we have two spectra 
for the bulge in the blue and red spectral range, respectively. 
However, the bars were observed only in the blue spectral 
range. 

Again, we have chosen two methods to analyse the 
stellar populations: the line-strength indices and a full 
spectrum-fitting approach. Unlike using the Lick/IDS sys- 
tem, we have chosen to perform our line-strength analysis on 
th e system imposed by t he ki nematics of th e galaxy, as done 
in lVazdekis et al.1 l|200ll ') and lYamada et al., (,200& ). For the 



second ap proach, w e use t he cr oss-correlation technique em - 
ployed bv lvazdekisi ()l999l l and lVazdekis fc Arimotol (jl999t ). 

In both cases, comparisons with stellar population syn- 
thesis models are needed; as for the kinema t ics, w e have 
used the empirical models of IVazdekis et"al] (l2010l). built 



from the MILES stellar library (|Sanchez-Blazauez et al.l 
l200d ). These models have a wide wavelength coverage (from 
3525 to 7500 A) and a resolution of 2.5A(FWHM) con- 
stant over the whole spe ctral range (jBeifiori et al.ll201ll : 
iFalcon-Barroso et al.ll201ll) . 



4-3.1 Emission line correction 

The ionized gas present in the galaxy leads to emission lines 
that contaminate the spectrum and fill in the absorption 
lines. Therefore, it is crucial to correct for this emission to 
get the properties of the stellar populations from the spectra. 
Indeed, if a proper separation of both contributions (the stel- 
lar absorptions and gaseous emissions) is not performed be- 
fore the analysis, the measured ages will be artificially older. 
In our particular case, we need to clean up the three Balmer 
lines included in our spectral range (H/?, H7 and H5), and 
the Fe5015 and Mgb indices, which can be affected by the 
[OIII]4959,5007 and [NII]5198,5200 lines, respectively. 

The common procedure for this correction is to mask 
the spectral region possibly contaminated by emission lines 
and fit the obtained spectrum with a set of stellar popula- 
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tion models. For th is purpose we use the ULyS^ package 
jKoleva et al.ll2009l ). which minimizes the between the 
observed spectrum and a model spectrum consisting of a 
linear combination of non-linear parameters (age and metal- 
licity), corrected for the kinematics and multiplied by a poly- 
nomial to avoid flux calibration problems. We subtract the 
ULySS best fit spectrum to the original one in order to get 
the residuals plus the contribution of the emission lines. We 
then approximate each emission line with a Gaussian, con- 
sidering only those features with an amplitude larger than 
three times the standard deviation of the residuals. 

Given the sensitivity of the emission correction, we 
double-check it by repeating it with the GANDALF package 
l|Sarzi et al.l 12009 ). The advantage of GANDALF is that it 
fits simultaneously the stellar and the gaseous contributions, 
by including emission lines as Gaussians in the set of stellar 
population models used as templates. Although GANDALF 
is also based on minimization, it differs from ULySS in 
that it fits the stellar contribution with a linear combination 
of models, while ULySS interpolates the grid of templates 
to look for the most appropriate single stellar population 
model. 

There are many details to take care of when fitting the 
full spectrum, such as the initial guesses and the number 
and range of the masked regions. For example, in the case 
of the blue spectra of NGC 357, we are forced to mask the 
CN band as we find evidences of its g-enhancement (no t 
included in the scaled-solar models of IVazdekis et~alll2010h . 
When the amount of emission derived from the two routines 
is in agreement, we subtract it from the observed spectra 
and use the results for the analysis of the stellar populations. 
Figure |4] presents an example of the fitting performed with 
ULySS for the bulge red spectrum because it shows H/J and 
[Olll]5007, which are the most prominent emission lines we 
observed. 

4-3.2 Line- strength indices 

A set of suitable absorption line-strength indices are mea- 
sured over the emission-cleaned spec tra of the bulge , inne r 
bar, and outer bar, as well as over the lVazdekis et all (|2010l ') 
models. In order to compare the results of the three compo- 
nents, the galaxy spectra and the SSPs are first degraded to 
their maximum resolution, given by the higher value of the 
velocity dispersion profile {a ~180 kms"'^). All the measure- 
ments are performed in the same system, i.e., the one defined 
by the kinematics of the galaxy in a consistent way. We plot 
the values obtained for an age indicator versus those of a 
metallicity indicator for the SSP models, in order to obtain 
model grids such as those shown in Figures [5] and |6l These 
grids are not perfectly orthogonal due to the age-metallicity 
degeneracy. The mean luminosity-weighted ages and metal- 
licities of each structural component of NGC 357 can then 
be derived by simply overplotting the observed indices on 
these grids. 

In particular, we computed all the Lick indices included 
in our spect r al ran ge, following the definitions given by 
iTraeer et al.l l)l998h . Note that although the index defini- 
tions are similar, our model measurements are not in the 

^ Available at http://uIyss.umv-Iyonl.fr/| 



Lick/IDS resolution and IDS instrumental system, but at 
constant a and fiux-calibration response curve. Since the red 
pseudo-continuum of the Fe4383 index falls partially out of 
our blue spectral range, we have defined a new Fe4383^^ 
index (named after Short Red, see Appendix]^ for details). 
Table [2] presents the definitions of the new Fe4383^^ and 
the original Lick/IDS Fe4383 index. 

For the bulge red spectrum, we compute the com- 
bined Iron-index (Fe) and total metallicity indicator [MgFejf] 
()Gonzalezlll993h . We avoid the use of the Fe5015 index, since 
it is highly contaminated by the prominent [OIII]5007 emis- 
sion line. 

In addition to the Lick age index-definitions for the 
Balmer lines (H/3, H7 and H^), we also me asure the 
more pow erful H/3o teervan tes fc Vazdekisll2009l ) and H7CT 
( Vazdckis fc Arimot o 1999) indices, which show great abil- 
ities to disentangle ages. The use of indices with enhanced 
sensitivity to the age or metallicity represents an advantage 
to lift the degeneracy between these parameters. 

4-3.3 Cross-correlation analysis 

The alternative technique we use to analyse the stellar popu- 
lations of NGC 357 is to cross-correla te the galaxy sp e ctrum 
with each SSP spectrum from the IVazdekis et al.l (|201G| ) 
model library. For this purpose the model spectra are pre- 
viously prepared to match the spectral range, velocity dis- 
persion, and spectral resolution of the data. Moreover, the 
galaxy and model spectra are rebinned logarithmically and 
normalized to remove the continua. 

In order to optimize the cross-correlation method for 
disentangling different stellar populations, it is necessary to 
adequately filter t he spectra and m ultiply them by a cosine- 
bell-like function ()Tonrv fc Davisl[i979i 'l . The importance of 
choosing a suitable filter lies in the possibility of getting rid 
of the noise in the spectrum; this purpose can be achieved 
by simply removing the largest wavenumbers, where the in- 
formation about the shortest wavelength ranges is included. 
Therefore, the limit is imposed by the resolution of the data. 
On the other hand, shorter wavenumbers contain informa- 
tion about wider spectral ranges, so possible residuals of 
the continuum removal due to errors in the fiux calibra- 
tion might also be filtered. The drawback of this procedure 
is that it implies a power loss of the final cross-correlation 
function, specially when filtering short wavenumbers where 
most of the signal is included. 

Apart from the filtering, it might be required to mask 
some regions in the original spectra, as is usually done in 
the full-spectrum fitting technique in the wavelength space. 
We tested different masks trying to avoid those features that 
are not well reproduced by the models due to mismatched 
abundance ratios: the CN in the blue spectral range and the 
Mg and H/3 features in the red spectral range. Again, these 
features contain most of the signal of the power spectrum 
so the choice of the masks has to be done very carefully in 
order to not lose most of the information. 

Finally, the peak height obtained for each correlation 
function is plotted against the model age and metallicity. 
Since the cross-correlation profile reaches a higher value 

3 (Fe) = (Fe5270 + Fe5335)/2, and [MgFe] = y'Mgb x (Fo) 
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Figure 4. The correction of tlie emission lines in tlie red spectrum of tlie bulge of NGC357 by means of ULySS. The observed and 
model spectra are plotted in black and orange, respectively in the upper panel. The vertical dashed lines mark the spectral regions we 
masked because of emission-line contamination (H/3 and [OIII]4959,5007 doublet) and a-enhancement (Mg). The residuals of the fit are 
shown in green in the lower panel. The Gaussians fitting the emission lines are plotted in magenta. 



Table 2. The definitions of the new Fe4383'5« and Lick/IDS Fc4383 indices. 



Index 


Blue pseudo-continuum 


Main bandpass 


Red pseudo-continuum 


Source 




(A) 


(A) 


(A) 




Fe4383'S'^ 


4359.125-4370.375 


4369.125-4398.087 


4419.328-4432.389 


this paper 


Fe4383 


4359.125-4370.375 


4369.125-4420.375 


4442.875-4455.375 


Traeer et al. (1998) 



when object and template are more similar (getting a maxi- 
mum value of 1 if they are exactly the same and no filtering is 
applied), the maximum value in the final peak heights plot 
indicates which template best resembles the galaxy spec- 
trum. This method takes advantage of the full information 
contained in the spectrum, instead of being constrained to 
certain features as in the case of the line-strength indices. 
This type of full-spectrum fitting technique works in the 
Fourier space and, as an additional advantage, it does not 
need for a previous emission correction of the data. The ap- 
plication of this technique to the der ivation of the st ellar 
populations was already int roduced bv lVazdekisI j 19991) and 
IVazdekis fc Arimotol l| 19991 ) and in this work it is adopted 
to double-check the results obtained from the analysis of the 
line-strength indices. 

To compute the errors for the cross-correlation method 
we perform a set of 100 Monte Carlo simulations by ran- 
domly adding Gaussian noise to the data, keeping the orig- 
inal S/N ratio. The 1-a regions are then overplotted to the 



cross-correlation peak height curves. In order to test the 
req uirements of th i s met hod we run some simulations with 
the lVazdekis et al.l l|2010l ) models by artificially changing the 
different parameters involved: spectral resolution, spectral 
range, and S/N ratio, taking as a reference the values mea- 
sured for our observed spectra. We find that the S/N of the 
data obtained for the outer bar is not enough to disentangle 
its age and metallicity, since the error bars for this case are 
too large (see Appendix [Bj . This prevents us from provid- 
ing a more constrained solution for the outer bar through 
the cross-correlation analysis. We have to note that the re- 
sults obtained here are aimed to be complementary to those 
from the line-strength analysis. Therefore, we restrict the 
application of the cross-correlation to the bulge and inner 
bar. These results are shown in Figures [7] and [8l and will be 
analysed in Sections 14.3.41 and 14.3.51 respectively. 
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Figure 5. Age indicators H/3o, H/3, H7ct=200i and }lSp (from top to bottom) versus the m etallicity indicators [M gFe] , Mgb, (Fe), CN2, and 
Ca4227, and the G-band (from left to right). The grids correspond to the SSP models of lVazdekis et al] 1I2OIOI) once smoothed to match 
the resolution of the data {cr ~180 kms~^). The lines represent different ages increasing from top to bottom (2.5, 3.2, 5.6, 10, and 18 Gyr, 
respectively) and metallicities, which increase from left to right ([Z/H] = —0.7, —0.4, 0.0, and 0.2, respectively). Ages and metallicities 
for the grids are indicated in the first panel. The black filled circle is the measure ment for the bulge of NGC 357, whereas the grey 
filled circles are the mean values for the bulges of a subsample of barred galaxies fr omlPerez e t al. (200(|) and lPerez &: Sanchez-Blazguej 
(|201lh . The open circles correspond to two elliptical galaxies in Virgo cluster from lYamada et al.l 1)2009 ) . All the reference galaxies have 
been selected to have a central velocity dispersion similar to that of NGC 357 [a ~180 kms"'^). 
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Figure 6. Age indicators H7o.=200 a^nd H^f (top and bottom panels, respectively) versus the metallicity indicators Fe4383'^^, CN2, and 
Ca4227, and the G-band (from left to right). The model grids are the same as in Figure [5] The coloured symbols are the measurements 
for the different structural components of NGC357: bulge (black circle), inner bar (orange diamond), and outer bar (blue triangle). The 
grey triangles are the mean values for the main bars of a subsample of galaxies from lPerez et al selected to have a similar central 

velocity dispersion than NGC 357 (cr ~180 kms~^). 



BULGE - RED SPECTRUM BULGE - BLUE SPECTRUM 
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Figure 7. Cross-correlation peak heights from the analysis of the red (left panel) and blue (right panel) spectra of the bulge of NGC 357 
with a set of SSP models of lVazdekis et al once smoothed to match the spectral range, resolution, and dispersion of the data. 

The different coloured lines correspond to metallicities of [Z/H] = 0.2 (black), 0.0 (red), —0.4 (blue), and —0.7 (green) as a function of 
age. The width of the coloured region represents the l-c uncertainty derived from Monte Carlo simulations. The largest cross-correlation 
peak height value in each plot is marked by an asterisk and the corresponding age with a vertical dashed line. 



4-3.4 Stellar populations of the bulge 

Figure [5]shows the analysis of the line-strength indices of the 
bulge of NGC 357. Four age indicators are plotted against 
several well selected metallicity indicators and the G-band, 
using features from both the red and blue spectral ranges 
available for this component. Balmer lines might be affected 
by emission line contamination, although the amplitude of 
this emission decays from the top (H/3o and H/3) to the bot- 



tom (H(5f) indices of Figure [S] In fact, we find no emission 
in the H5 feature for the bulge. For the case of H/3o and H/3, 
where the emission contamination is larger, we find some 
disagreements in the resulting ages with respect to those 
derived from H7 and H(5. Although the emission correction 
has been performed very carefully (see Section [4.3.11) . these 
inconsistencies are probably due to an underestimation of 
the emission affecting the H/3 feature within the noise level. 
As an attempt to take this effect into account, the H/3 and 
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Figure 8. Same as in Figur^Jfor tlie blue spectrum of the inner bar of NGC 357. 



H/3o error bars are obtained by computing the maximum 
and minimum (but no null) 1-a emission correction. 

Moreover, we find subtle differences between the ages 
obtained with the H/?o and H/3 indices. These differences 
might be attributed to either the emission correction, that 
affects the two definitions in a different manner, or to the 
greater sensitivity of H/3o to the a-enhancement in com- 
parison to that of H/3 (|Cervantes fc Vazdekisll2009l '). Indeed, 
these differences are more evident in the a-enhanced indi- 
cators (Mgb and CN2), whereas the results from (Fe) or 
Ca4227 are almost identical. Taking into account the error 
bars and the results from all the plots in Figure [H we derive 
a mean luminosity-weighted age of ~8 Gyr for the bulge of 
NGC 357. 

From the first column in Figure [SI where the total 
metallicity indicator [MgFe] is plotted, we derive a super- 
solar metallicity ~0.2 dex for the bulge of NGC 357. For 
this component we also measure the (Fe), Mgb, and CN2 
indices, so that it is possible to estimate the [Zjvig/Z^pc)] 
and [ZcN2/Z(Fe>] ratios, which are calculated as the rela- 
tive metallicity difference obtained when plotting the various 
metallicity indicators versus the Mg (|Vazdekis et aLll200j ). 
[ZMg/Z(Fe)] and [ZcN2/Z(Fe>] are shown to be good prox- 
ies for the [Mg/Fe] and [CN/Fe] abundance ratios. In fact, 
fZ Mgb/Z/Fe\] and [ Mg/F e] follow a linear relation, as shown 
in lVazdekis et all (|2010l '). These overabundances are impor- 
tant because they can be used to constrain the timescale for 
the formation of the bulk of their stellar populations: Mg 
(as well as other a elements) is the product of type II super- 
novae, which explode soon after the star formation; on the 
contrary, the Iron-peak elements are ejected by type la su- 
pernovae, on longer timescales of around I Gyr. In the mean- 
time between those two cases, the low- and intermediate- 
mass stars enrich the interstellar medium with C and N. 
Therefore, the [ZMgt/Z^Fc)] and [Zcn2/Z(fc>] values might 
help to put constraints on the assembly history of a galaxy. 
In the case of NGC 357, we find [ZMgb/Z(Fo>]=0.3-0.4 dex, 
which is the value expected for an E/SO galaxy of a ~180 
kms~^ according to the well-known correlation between 
[ZMg/Z^Fc)] and a for early-type objects (e.g., iJ0rgenseni 



199£ ; iTrager et al.l I2OO0I: iKuntschneij I2OO0I : IVazdekis et al] 
2004 : Icarretero et al.ll2004l '). Note that an extrapolation of 



the grid is necessary to estimate the [Mg/H] abundance, 
since the models do not extend beyond [Z/H]=0.2 dex. We 
also find a supersolar [Zcn2/Z(fc)] ratio for the bulge of 
NGC 357, that reaches ~0.2 dex. 

Note that the bulge of NGC 357 shows subsolar values 
for the Ca4227 index. This is the expected behaviour for 
galaxies of similar central velocity dispersion. In fact, rather 
than resembling other a ele ments as th e Mg, Ca tends to 
get subsolar values tVazdekis et al.lll997f ). 

In order to find out if the formation of the bulge of 
NGC 357 has followed a distinguished process due to the 
presence of the inner bar, we compare our results with 
those obtained for various elliptical galaxies and other spi- 
ral bulges. Thus, we have overplotted the index values for 
two Virgo elliptical s with a centra l a ~ 180 kms^^, taken 
from the sample of lYamada et al] (|2006h . These measure- 
ments appear as two open circles in Figure [5] Moreover, 
we have selected a subsam ple of four early-type b arred 
galaxies from the works by iPerez et all (|2007l . |2009| ) and 
iPerez fc Sanchez-BlazQue3 (|201ll ) with central velocity dis- 
persions similar to that of NGC 357. Two galaxies host a 
single bar (NGC 1169 and NGC 1358), and two are double- 
barred systems (NGC 2859 and NGC 2962). For these ob- 
jects, we have measured the mean line-strength indices in- 
side the bulge region from the values at a resolution of 
a ~180 kms~^. The results, when possible, are also shown 
in Figure [5] as grey filled circles. In general, the age, metal- 
licity and a-enhancement obtained for NGC 357 resemble 
very closey those of the comparison bulges. The age and 
metallicity of the Virgo ellipticals are also in agreement with 
NGC 357, although they show a larger [Zcn2/Z(fc>] enhance- 
ment, with a mean value above ~0.3 dex. 

The results for NGC 357 from the line-strength analy- 
sis are also verified by the cross-correlation method. Figure 
[7] shows the peak heights resulting from cross-correlating 
the bulge spectra with SSPs of different ages and metallic- 
ities, following the methodology described in Section [4.3.31 
From both spectral ranges, we obtain a supersolar value 
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([Z/H]'^0.2 dex) for the mean luminosity- weighted metal- 
licity, and a final age of 7-8 Gyr (although the error regions 
span a wider range). So far, the cross-correlation method 
alone is not able to clearly distinguish a mean luminosity- 
weighted stellar population above ~5 Gyr, due to the flat- 
tening of the cross-correlation peak heights versus age pro- 
files and the amplitude of the error bars. For younger ages. 
Figure [7| shows that the curves are steeper, allowing to dis- 
entangle between different populations very efficiently, ffow- 
ever, the fact that the results obtained from this proce- 
dure directly agree with the line-strength analysis, points 
out the potential power of this method. Note that the age- 
metallicity degeneracy is evident in Figure [T] the shoulder 
of each metallicity curve (i.e., the point where there is a 
clear transition in the slope) occurs at older ages for lower 
metallicities. 



4-3.5 Stellar populations of the inner and outer bars 

The relative ages and metallicities between bulge, inner bar 
and outer bar of NGC 357 are shown in Figure |6l where 
the age indicators H7o- and H5f are plotted against various 
metallicity indicators (including the new Fe4383^'^ defini- 
tion) and the G-band. In this case we are limited to indices 
in the blue range, since it is the only spectral range available 
for the two bars. 

In the first panel of Figure [S] we plot Il7o- against 
Fe4383^^, which provide the most orthogonal grid thanks 
to the age-disentangling power of 117^. We find that the 
bulge and inner bar show approximately the same mean 
luminosity-weighted age (~8 Gyr) and solar metallicity. The 
outer bar also presents a similar age as the inner structures, 
but it tends to be less metal-rich ([Fe/H]'^-0.4 dex). How- 
ever, the error bars are large due to the S/N requirements of 
the H7ct index. For this reason it is necessary to complement 
these results with those from the less orthogonal HSf ver- 
sus Fe4383^^ grid, which confirms the previous estimations 
for the bulge and inner bar. The trend towards lower metal- 
licites for the outer bar is also verified. In fact, the only 
difference we find between the two plots is that the outer 
bar seems to be slightly older when looking at the H5f ver- 
sus Fe4383^^ grid. However, taking into account the error 
bars and the results from the rest of the plots (particularly 
USf versus CN2 and H7CT versus Ca4227), we must conclude 
that the three structural components of NGC 357 are nearly 
coeval. 

Since we have measurements for the Fe4383'^^ and GN2 
indices, we can use the [Zcn2 /Zpo43g3SR] ratio as a proxy 
for the [GN/Fe] abundance, as we did in Section [4.3.41 with 
[ZcN2/Z(Fo>]- We find that, whereas the bulge and inner bar 
present an overabundance of [ZcN2/Zpe43g3SR]>0.2 dex, the 
outer bar has a more enhanced value above 0.4 dex. 

To make a comparison with the age and metallicity of 
other main bars, we take into account the same galaxies we 
adopted as a reference in Figure [5] We compute the mean 
values of the line-strength indices along the m ain bar at a 
resolu tion of a ~f80 kms^^ from the data of iPerez et al.l 
These comparison values are shown in Figure [6] and 
the plotted errors take into account the variation of the line- 
strength indices along the main bar. Like for the bulge case, 
the outer bar of NGG 357 shows no particular features when 



compared to the other main bars. Indeed, its age and metal- 
licity are within the ranges defined by the control sample. 

We apply the cross-correlation method only to the in- 
ner bar spectrum since that of the outer bar has too low 
S/N. The mean luminosity- weighted age and metallicity are 
consistent with those derived from the analysis of the line- 
strength indices. Figure[8]shows that although the maximum 
of the cross-correlation peak heights is obtained for ~4.5 
Gyr, the solution spans over a wider range (from about 4 to 
10 Gyr) when taking into account the error regions. The re- 
sulting metallicity, despite the partial overlapping between 
the solar and supersolar curves, is [Z/H]~0.2 dex. 



5 DISCUSSION 

The various pieces of evidence provided by this work lead to 
two main questions: which is the real structure of NGG 357, 
and how this galaxy has been formed. The analysis of the 
photometry and kinematic constrains the structural compo- 
sition of NGG 357, which reveals itself as a really complex 
galaxy. On the other hand, the study of the stellar popu- 
lations provides surprisingly homogeneus results and sheds 
light on the most probable assembly process for our galaxy. 
For the sake of clarity, we split the Discussion in two Sec- 
tions, devoted to each of these questions. 



5.1 Structural components of NGC 357 

Photometrical studies have confirmed that NGG 357 is com- 
posed of, at least, fo ur structure s , namely bulge, inne r bar, 
outer bar and disc (|Erwinll2004 Uguerri et al.ll2005l ). The 
presence of the inner bar is also supported by the photo- 
metrical and kinematical analysis performed in this work. 
In fact, the ellipticity and PA profiles presented in Figure 
[2I and the velocity and velocity dispersion profiles along the 
inner bar direction presented in Figure [3] show clear signa- 
tures of this small, secondary bar, as explained in Sections 
[4TI and 14:21 

Interestingly, the analysis of the kinematics also shows 
evidence of an additional, inner component: a kinematically 
decoupled structure appears at the centre of NGG 357, ro- 
tating faster than its surroundings. This decoupling spans 
approximately ±2 arcsec and matches in size with the ob- 
served central cr-drop. The different structures that shape 
NGG 357 have to account for all the signatures present in 
the kinematical profiles. In this Section we state the two 
possibilities that match with the results obtained in this 
work. 



5.1.1 An inner disc and a classical bulge 

The presence of a central velocity dispersion minimum in 
a spiral galaxy was first reported bv lBottemal (|l989l 'l. Since 
then, the number of g alaxies showing a g-drop has increased 
significantly. In fact, iFalcon-Barroso et al.l (|200m find this 
signature in at least 46% of a sample of 48 early-type spi- 
rals, alth ough other cases have been noticed in smaller sam- 
ples fe.g. jMarguez et al.ll200j ). particul arly in barred galax - 
ies and even including late -type spirals (Ganda et al.ll200g 'l 
and double-barred objects (|Emsellem et al.ll200ll ). The most 
accepted explanation for the formation of a"-drops in disc 
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galaxies is the star formation at their central regions: the 
new stars acquire the kinematical properties of the gas 
they are formed from; the dissipative nature of that gas is 
more efficient in the central regions because of the higher 
density and converts the gas into a cold stellar compo- 
nent, such as an inner disc, with a lower velocity disper- 
sion than the surroundings. This scenari o is in agree ment 
with the N-body simula t ions o f .Wozniak et alj l|2003l ) and 
IWozniak fc ChampavertI (|2006l ). who made use of bars to 
transport gas to the central regions. 

The case of NGC 357 seems to support that hypothesis 
since the central decoupling in velocity might be due to an 
inner disc whose low velocity dispersion would be the cause 
of the (j-drop. Therefore, NGC 357 is composed of at least 
five structures: bulge, inner disc, inner bar, outer bar and 
disc. 

The CT-drop means a local minimum in the velocity 
dispersion with respect to the corresponding value for the 
bulge, which reaches a ~180 kms^^ or higher (the velocity 
dispersion profile should peak at the centre, where the a- 
drop masks the velocity dispersion of the bulge). This high 
velocity dispersion value for the bulge implies it is pressure 
su pported and cl assical. This result is supported by the work 
oflAguerrTelaD ([2005), who find that the bulge of NGC 357 
follows the same fundamental plane as the ellipticals or clas- 
sical bulges of other early-type objects. 

Within this context, the appearance of the cr-hoUows 
at the edges of the inner bar in the velocity dispersion pro- 
file is also well understood. The cr-hollows were first seen 
in a sample of four double-barre d galaxies analysed throug h 
integral-field spectroscopy ()de Lorenzo-Caceres et a"l]|2008l ). 
By means of N-body simulations, they tested the possible 
scenarios that might give rise to these hollows, showing that 
their size and amplitude depend on two parameters: the 
relative contribution of the flux of the bulge to the total 
luminosity and the difference in the velocity dispersion of 
the two components. Therefore, a galaxy hosting a classical 
bulge with a high velocity dispersion and a colder, inner bar 
which dominates the total luminosity at its ends (where the 
light profile of the bulge has already decayed) will clearly 
show the two cr-hollows at the bar edges, as for the case of 
NGC 357. In this scenario, the bulge is the hottest structural 
component of NGC 357 and the velocity dispersion decreases 
for the outer structures. Table[3]indicate the velocity disper- 
sion values for each component. It is important to note that 
the actual velocity dispersion of the inner disc is probably 
lower than the value given in Table [3l since the measured 
value includes the contribution of the bulge within the inner 
disc region. 

There are two important drawbacks for this possible 
structural composition of NGC 357. First, the a4 > signa- 
ture of the inner disc should be detected in the photometri- 
cal analysis presented in Section 14.11 This is the case of the 
nuclear stellar discs observ ed in the centre of ellipticals and 
bulges (|Morelli et al.|[2O10l . and references therein) . A careful 
inspection to the inner ±2 arcsec in Figure [2] shows no evi- 
dence of discy isophotes. Second, the formation of the inner 
disc is usually related to a recent star formation in the cen- 
tre, but the analysis of the stellar populations for NGC 357 
shows old ages for the central components; a significant 
young population such as that of an inner disc should be 
noticed when studying the mean luminosity-weighted age. 



even if it is mixed with an older, more massive component. 
Therefore, although this hypothesis seems promising, it does 
not explain all the properties found in this work. 

5.1.2 A pseudobulge 

The second possibility for the structure of NGC 357 is that 
there is no inner disc. If this is the case, the central kine- 
matical decoupling and the a-drop have to be caused by the 
bulge itself, which therefore shows disc-like properties that 
indicate it is a pseudobulge rather than a classical on e. This 
hypo thesis is strongly supported by the Sersi c index (ISersic 
'1968") measured for the bulge of NGC 357 by lAguerri et al. 
(2005 ); they performed a photometric decomposition of an 
/-band image and obtained n=1.4 0±0.08, which i s com - 
patible with a pseudobulge. Indeed, iFisher fc Drorvl (|2008l ) 
study a sample of 77 galaxies and find that ~90% of pseu- 
dobulges have n<2, whereas classical bulges have n>2. 

Within this new context, the cr-drop is not actually a 
drop with respect to the higher velocity dispersion of the 
bulge, but the maximum velocity dispersion of the pseu- 
dobulge itself. Of course, the pseudobulge presents a not so 
high velocity dispersion, since it is rotationally supported, 
and therefore it can be colder than the inner bar. In this 
case, the two peaks of the velocity dispersion values at ei- 
ther sides of the centre have to be understood as signatures 
of the relatively hotter inner bar where this component dom- 
inates the total luminosity. The cr-hollows appear where the 
bar ends and the light of the pseudobulge is the main contri- 
bution to the total flux. The surrounding structure, i.e., the 
outer bar, is also hotter than the outer parts of the pseu- 
dobulge and for that reason the velocity dispersion values 
increase again beyond the cr-hollows. The characteristic ve- 
locity dispersions for each structure is indicated in Table |3l 
This scenari o is completely different from tha t of the galax- 
ies shown in lde Lorenzo-Caceres et al. 1(2008). However, the 
main conclusion on the origin of the cr-hollows remains as it 
is attributed to the contrast of velocity dispersions between 
components. 

This hypothesis also explains the results from the anal- 
ysis of the morphology of NGC 357: the ellipticity profile 
acquires approximately the same value at either sides of 
the peak due to the inner bar, suggesting that the centre 
and the surrounding regions are dominated by the same 
structure. Concerning the dynamical support of the bulge, 
we estimate Vmax/cr ^^0.23 and e<0.1 in the bulge region; 
therefore, the bulge i s consistent with an oblate sp heroid 
fiattened by rotation (|Kormendv fc Kennicuttl |2004| ) . Nev- 
ertheless, the interpretration of the Vmax/cr — e relation has 
to be done carefully, since a complex kinematical structure 
of the ga laxy might conduct to misleading results, as it is 
shown in lEmsellem et all l|201lh . 

5.2 Formation of NGC 357 

The analysis of the stellar populations of NGC 357 is focused 
on three single regions: bulge, inner bar and outer bar. The 
striking result obtained in this work is that the bulge and in- 
ner bar show very similar stellar population properties, indi- 
cating that they were formed in the same process or, at least, 
from the same source of stars. Moreover, the outer bar is also 
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Table 3. Characteristic velocity dispersion values for the different structural components of NGC357 within the context of the two 
possible scenarios. For the bulge and pseudobulge, these values represent their maximum velocity dispersion. For the outer bar and disc, 
the mean velocity dispersion values within the corresponding regions included in Figure[3]are provided. 



Scenario 


(Pseudo)bulge 


Inner disc 


Inner bar 


Outer bar 


Disc 


(1) 


>180 


~170 


~165 


~145 


~105 


(kms-i) (2) 


-170 




~180 




~105 



(1) Classical bulge + inner disc + inner bar + outer bar + disc; 

(2) Pseudobulge + inner bar + outer bar + disc. 



nearly coeval to the inner structures, although it presents a 
lower metallicity and a larger a-enhancement. This result 
points towards a faster assembly of the outer bar with re- 
spect to the bulge and inner bar. Finally, we compare the 
bulge and outer bar of NGC 357 with other bulges and main 
bars of single- and double- barred galaxies with similar cen- 
tral c r values (iPerez et al ] 12009; Porcz fc Sanchez-Blazquej 
I2OIII ). finding no significant differences among them. 

The results obtained here for NGC 357 indicate that 
theories claiming that inner bars play a major role in the 
secular evolution might not be suitable for this galaxy. If in- 
ner bars were key to transport gas to the central regions and 
trigger star formation there, a gradient in age and metallicity 
from outer to inner structures would be expected, present- 
ing a younger population in the galaxy centre. Although our 
galaxy has gas content and probably some star formation is 
taking place in it, the stellar population analysis indicates 
that this is not the main mechanism driving the evolution 
of NGC 357. Therefore, these results lead to the question on 
how NGC 357 formed. 

Numerical simulations have predicted different sce- 
narios for the formation of nested bars systems, most 
of them requiring the presence of a dissipative compo- 
nent that is dy namically disturbed and generates the sec- 
ondary bar (e.g..lFriedh fc M artinet" 1993': 'Hell er et al.ll200ll : 



IShic 



Helled l2003 ~iEnglmaicr & Shlo smanl |2004| ). 



Within this scenario, the inner bar should be younger than 
the outer bar. However, the two bars of NGC 357 have 
similar population properties, so they were formed at the 
same time or, most likely, the existing stars from the outer 
bar or the disc redistributed and shaped the inner struc- 
tures. Therefore, only those numerical works that do not 
need gas to generate the double-barred objects seem to be 
compatible with the observational evidences shown here for 
NGC 357. This is t he ca se of the N-body simulations by 
iDebattista fc ShenI |20o3), who use a disc and a rapidly ro- 
tating bulge (like a pseudobulge) to create a double-barred 
system. In these simulations, the inner bar appears first, but 
this is probably a consequence of the initial presence of the 
pseudobulge: in the real Universe, the outer bar most likely 
is formed first, so the pseudobulge would be generated from 
the gas flow along it. This scenario reinforces the hypothesis 
of the pseudobulge for the case of NGC 357, since the bulge 
rotation is an essential condition for generating the two bars 
within the IDebattista fc ShenI (|2007l ') gas-free framework. 

There is one possibility to reconcile the classical sec- 
ular approach with the observational results for NGC 357, 
and this is to assume a very fast formation of the central 
components: immediately after the outer bar was created, 



the gas flowed along it, triggered a quick star formation and 
formed the inner bar and the bulge. Within this scenario, 
the flnal ages of the three components would be very simi- 
lar. The outer bar would be slightly older, but this is con- 
sistent with our results taking into account the error bars. 
In this scenario a positive metallicity gradient for the gas 
along the outer bar is required. Although this fact might 
be somewhat striking, from a stellar (not gaseo us) point of 
view there is evidence for such gradients. Indeed. TFerez et al.l 
(2009) study a sample of 20 barred galaxies and flnd that 
some of them tend to become more metal-rich at the ends of 
the main bars. Moreover, those objects have older mean ages 
and higher central velocity dispersions than the rest of the 
sample, which is in agreement with the properties found for 
NGC 357. Unfortunately, the S/N ratio of our data is not 
enough to measure gradients along the main bar to check 
this possibility. 



6 CONCLUSIONS 

We have performed for the flrst time a detailed analysis of 
the morphological, kinematical, and stellar population prop- 
erties of a prototype double-barred early-type disc galaxy. 
We put special emphasis on the inner bar in order to con- 
strain its role in the evolution of this galaxy. The observa- 
tional strategy consisted in taking very deep long-slit spec- 
tra along the directions of the inner and outer bars and the 
major axis of the disc. 

The presence of the two bars is clearly revealed by 
the ellipticity and FA profiles, which show the usual 
sharp changes due to this kind of structure. Moreover, 
the kinematical analysis also presents clear signatures due 
to the inner bar: the double-hump profile in the velocity 
jEmsellem et al.ll2001^ and the g-hollow s in the velocity dis- 
persion ( de Lorenzo-Caceres et a"l]|2008l ). These hollows are 
due to the contrast between the velocity dispersion values 
of the inner bar and the central structure, that might be 
a classical bulge with a velocity dispersion higher than the 
inner bar, or a pseudobulge with a lower velocity dispersion 
than the bar. The relative contribution of the inner bar to 
the total fiux also influences the size and depth of the a- 
hoUows. It remains unclear if the cr-hoUows or other related 
features may appear for the case of the outer bar. Unfortu- 
nately, the S/N of the main bar spectrum is not sufficient 
to explore the kinematics at its ends. 

The analysis of the kinematics of NGC 357 also reveals 
a kinematically-decoupled structure rotating faster than its 
surroundings at the centre of the galaxy. This signature 
matches in size with a central tj-drop. NGC 357 is the first 
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observed galaxy in which a cr-drop and the tj-hollows coexist. 
This fact illustrates the differences between both signatures: 
whereas the former are just due to the presence of the in- 
ner bar, the latter is caused by an inner, colder structure. 
Due to the different nature of both signatures, we want to 
stress that a double-barred galaxy might present a a-drop, 
the (T-hoUows, or both, as in this case. 

The structural composition of NGC 357 has to account 
for all the observational evidences found in the analysis of 
the photometry and kinematics. The presence of the disc and 
the two bars is clear, but the central region turns out to be 
more difficult to disentangle. In fact, two possible scenarios 
are compatible with the results obtained in this work. The 
first possibility is that NGC 357 hosts a classical, hot bulge 
together with an inner disc, which would be responsible of 
the decoupling in the velocity profile and the cr-drop. The 
second scenario is that the bulge of NGC 357 is, in fact, a 
cold pseudobulge, which is directly related to the kinemati- 
cal central signatures. 

The main conclusion of this work is that the bulge 
and inner bar show similar stellar population properties 
(age, metallicity and a- enhancement), whereas the outer 
bar has no significant difference in age but it is less metal- 
rich and more a-enhanced, indicating that it was assembled 
in shorter timescales than the inner structures. This result 
seems to discard for this galaxy the traditional secular evo- 
lution scenario, in which the star formation triggered by 
the gas flown along the outer structures is causing the for- 
mation of the central parts. Therefore, NGC 357 has been 
shaped by the redistribution of the existing stars, maybe 
in a secular way in which stars from the disc form the in- 
ner components, or maybe in an initial disc + pseudobulge 
framework so the inner bar is lately formed from the pseu- 
dobulg e stars. This last scenario is backed by the numerical 
work of lDebattista fc ShenI fioO?) , who create double-barred 
systems from a pure stellar disc and pseudobulge. In this sce- 
nario, the fast rotation of the pseudobulge is the key element 
to generate bars without the gas contribution. 

An analysis such as the one presented in this work for a 
larger sample of double-barred galaxies is required to derive 
general and robust conclusions on the formation of these 
complex objects. Fortunately, the striking results obtained 
for NGC 357 are motivating and shed light on the impor- 
tance of carefully studying double-barred galaxies, from the 
outer disc to the very central regions. 
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APPENDIX A: THE NEWLY DEFINED 
FE4383^^ INDEX 

For the inner and outer bars of NGC 357 we only have blue 
spectra, from 3990 to 4440 A. There is only one Iron index in 
the Lick system within this range, at 4383 A. Unfortunately, 
the red pseudo-continuum of this Lick Fe4383 index falls par- 
tially out of our blue spectra, so we have redefined it in order 
to introduce a new Fe4383'^" index suitable for our spectral 
range requirements. Fe4383^^ is named after Short Red (due 
to its motivation) and keeps the blue pseudo-continuum of 
the Lick Fe4383 index, but modifies the range definitions 
for the main bandpass and the red pseudo-continuum (see 
Table [2]) . This new index has been tested with the bulge 
of NGC 357, for which we have data in two spectral ranges 
so we can measure additional Iron indices. Figure lATI shows 
several age indicators plotted against Fe4383'^^ and (Fe). 
The corresponding measurements for the bulge of NGC 357 
are overplotted. The age and metallicity obtained from the 
two Iron index definitions are consistent. We therefore use 
Fe4383'^^ for estimating the metallicity in our blue spectra. 



APPENDIX B: S/N CONSTRAINTS FOR THE 
CROSS-CORRELATION METHOD 

As explained in Section 2331 we explored the requirements 
of t he cross-correla t ion m ethod through several tests with 
the IVazdekis et~all (|2010l ') model library, by changing their 
spectral resolution, spectral range and S/N ratio. We found 
that the S/N of the outer bar spectrum is not sufficient to 
disentangle its age and metallicity. This fact is illustrated in 
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Figure Al. Comparison of tlie results obtained for tlie bulge 
of NGC 357 using as metallicity indicators the newly defined 
Fe4383'^^ and the (Fe) indices. The model grids are the same 
as in Figure [S] 



Figure [BTI where we show the results from aplying the cross- 
correlation method to a model of 8 Gyr and solar metallicity, 
once cut and smoothed to mimic the observed spectrum for 
the outer bar. We also added gaussian noise to the model to 
match the S/N ratio of the data. The error regions are so 
large that, for ages above ~3 Gyr the results for the solar 
and supersolar metallicity are completely mixed. 
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